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Abstract. Using optical/near-IR broadband photometry together with Ho- emission line data, we attempt to constrain the star 
formation histories, ages, total stellar masses and stellar mass-to-light ratios for a sample of extremely blue low surface bright- 
ness galaxies. We find that, under standard assumptions about the stellar initial mass function, the Ho- equivalent widths of 
these objects appear inconsistent with recently suggested scenarios including constant or increasing star formation rates over 
cosmological time scales. In a critical assessment of the prospects of obtaining ages from integrated broadband photometry, 
we conclude that even with near-IR data, the ages are poorly constrained and that current observations cannot rule out the 
possibility that these objects formed as recently as 1-2 Gyr ago. Methods which could potentially improve the age estimates 
are discussed. The stellar masses of these galaxies are inferred to lie below 10'° M G . This, in combination with low ages, could 
constitute a problem for current hierarchical models of galaxy formation, which predict objects of this mass to form predomi- 
nantly early in the history of the universe. The possibility to use the ages of the bluest low surface brightness galaxies as a test 
of such models is demonstrated. 
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1. Introduction 

Low surface brightness galaxies (LSBGs), here characterized 
by central surface brightness levels fainter than figfi ~ 23 
mag arcsec~ 2 , have stellar populations ranging from very red 
(B - V ~ 1.2 mag) to very blue (B — V ~ 0.2 mag) (e.g. 
O'Neil et al. 119971 1. The blue LSBGs are among the most 
gas-rich galaxies in the local universe (e.g. McGaugh & de 
Blok I1997t and have low metallicities (e.g. McGaugh [19941 
Ronnback and Bergvall 1 19951 1. typically 5-20% of the solar 
abundances. These properties could suggest that such galax- 
ies, of which most have disk-like structures, may have formed 
recently. 

According to currently popular models of structure forma- 
tion based on the cold dark matter (CDM) paradigm, low-mass 
galaxies should predominantly form early, whereas higher- 
mass objects are predicted to form through mergers and accre- 
tion at later epochs. This hierarchical scenario for galaxy for- 



mation has recently been questioned by the discovery of anti- 
hierarchical relations between mass and mean stellar age for 



low-redshift galaxies (e.g. Thomas et al. ;2002 Jimenez et al. 
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2005 ) and the discovery of large populations of high-mass ob- 
jects present already at high redshift (e.g. Somerville 2004 1. 
If the existence of truly young, low-mass galaxies could be 
demonstrated, this could turn out to be equally problematic. In 
this context, discussions on individual candidates for galaxies 
having formed recently, such as SBS 0335-052, Tololo 65 and 
IZwl8 (e.g. Izotov & Thuan fT999l Kunth & Ostlin 120001 are 
however of limited value, since even in the hierarchical picture, 
a few stragglers are expected (e.g. Nagamine 2001 ). Only the 
discovery of a substantial number of young, low-mass galaxies 
would impose a serious constraint. Since LSBGs are believed 
to contribute significantly to the galaxy number density of the 
universe (e.g. McGaugh et al. [191)51 O'Neil & Bothun l2TjQ0l 
and since many of these objects appear to be both unevolved 
and to have low masses, it becomes important to establish ex- 
actly when they formed. 

The ages of the LSBGs have long been a matter of con- 
troversy. In a search for young galaxies in the local universe, 
Ronnback and Bergvall ( 1994) carried out one of the first mul- 
ticolour studies of these objects. The sample used was se- 
lected on the basis of colour criteria, favouring extremely blue 
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galaxies with only small radial colour gradients. Ronnback and 
Bergvall (1995 1 derived chemical abundances for part of this 
sample and found very low metallicities. This was argued to be 
the major reason for the blue colours, as opposed to low ages 
- which was the most popular interpretation at the time. In a 
preliminary study (Bergvall & Ronnback [19941 . where BVI 
photometry was interpreted using spectral evolutionary mod- 
els, ages typically higher than 3 Gyr were found. Padoan et 
al. J1997I and Jimenez et al. (1998 1 similarly used UBVRI 
photometry to conclude that most of the galaxies in their sam- 
ple appeared to be older than about 7 Gyr. The galaxies stud- 
ied were however not extremely blue, and tendencies for ages 
below 7 Gyr for the bluest galaxies in their sample were in- 
deed present. Schombert et al. (2001 1 argued on the basis of 
their V - I colours and relative HI content that the most gas- 
rich LSBGs should typically have mean stellar ages below 5 
Gyr. Possible effects of gas infall, which could make the galax- 
ies appear younger than their true age, were however not dis- 
cussed. An opposing view was stated by van den Hoek et al. 
(2000 1 who, based on UBVRI data, concluded that the gas- 
rich LSBGs appeared to be very old disks with a strongly sup- 
pressed star formation activity, although the possibility that the 
bluest LSBGs might be dominated by stellar populations much 
younger than 5 Gyr could not be ruled out. 

A serious shortcoming of these investigations is their re- 
liance on optical photometry. Optical colours are, by them- 
selves, blunt tools for constraining the ages of blue LSBGs, 
since the youngest stars may dominate the optical part of the 
spectrum. An old stellar component will thereby be extremely 
difficult to detect, unless strong population gradients allow the 
selection of regions in the galaxies where the old, red popu- 
lation contributes significantly to the integrated optical light. 
With the inclusion of near-IR broadband photometry one would 
naively expect to be able to constrain the ages more strin- 
gently. An investigation of this type was carried out by Bell 
et al. (20001 for a sample of LSBGs covering a wide range 
of colours. Although a trend of lower mean stellar ages with 
bluer colours was found, a simplifying assumption adopted in 
this work was that the galaxies studied all formed 12 Gyr ago, 
making the average stellar age a function of metallicity and star 
formation history only. 

Since LSBGs are believed to be more dark matter dom- 
inated then their high surface brightness counterparts (e.g. 
McGaugh & de Blok fT998l Zavala et al. l2003l . they have be- 
come prime targets for testing predictions of dark halo density 
profiles through rotation curve decomposition (e.g. de Blok et 
al. I2001i . In this procedure, some mass-to-light ratio (M/L) 
believed to be representative of the stellar component in ques- 
tion is typically assumed. The most appropriate value for this 
parameter, which is intimately linked with both age and star 
formation history, is however still uncertain. It is therefore im- 
portant to determine the range of M/L allowed by the observed 
properties of the stellar populations of LSBGs. 

In this paper, optical/near-IR photometry and Ha data for 
a sample of 9 extremely blue LSBGs are analyzed to con- 
strain the properties of the stellar populations in these objects. 
Spectral evolutionary models covering a wide range of param- 



eter values at high resolution are used to investigate their star 
formation histories, ages, stellar M/L and total stellar masses. 

2. The data 

The optical/near-IR photometric data for the 9 blue LSBGs in- 
vestigated here have previously been presented in Ronnback & 
Bergvall dl994l and BergvaU et al. ( fT^99l . 

The targets used in the optical study of Ronnback & 
Bergvall ( 1994) were selected from the ESO-Uppsala cata- 
logue (Lauberts & Valentijn 1989) with the criteria of: 

1) A mean surface brightness inside the effective radius of 
/u e g > 23.5 mag arcsec" 2 ; 

2) An average colour <B - R> of less than 0.5 mag in re- 
gions with surface brightness //g = 20.5-26 mag arcsec -2 ; 

3) A(B - R) < 0.2 mag, where A(B - R) represents the 
difference in colour inside and outside the effective radius. 

In the near-IR photometric study by Bergvall et al. ( 1999 1, 
the surface brightness criterion was however slightly relaxed. 
Here, galaxies with a surface brightness 0.5 mag brighter than 
the original study, i.e. ff^ > 23.0 mag arcsec -2 , are therefore 
considered. 

The colour criterion targets the bluest 5% of the LSBGs 
(j/gf > 23.5 mag arcsec -2 ) in the ESO-Uppsala catalogue. The 
galaxies selected this way are dust-poor (Bergvall et al. 1999 1, 
which reduces the uncertainties in extinction corrections when 
using models to interpret their properties. The requirement of 
small colour gradients avoids situations where two or more sep- 
arate stellar populations are located in different regions of these 
objects, which further simplifies the analysis. 

The data are summarized in Table \l\ The oxygen abun- 
dances, O/H = 12 + log(/Y(0)/jV(H)) are from Ronnback 
and Bergvall ( 1995 1. The redshift and HI mass of 074-16 has 
been adopted from Kilborn et al. (2002i. In all other cases, the 
HI data stem from Bergvall et al. (2005). The Ha equivalent 
widths, EW(Ha), of 3 galaxies are derived from narrowband 
images obtained at the ESO/MPI 2.2 m telescope, La Silla, and 
at the NOT 2.5 m telescope on La Palma. The estimated lcr 
errors are cr(B - V) - 0.04 mag, cr(<B - R>) = 0.10 mag, 
cr(y - = 0.06 mag, cr(B - J) = 0.08 mag, cr( J - H) = 0.12 
mag, cr(H - K') = 0.14 mag, cr(0/H) = 0.1 and o-(EW(Ha)) = 
20%. 

In Table [fl both the central surface brightness of the ex- 
ponential disk and the "true" surface brightness of the centre 
are presented. This is particularly useful in the cases where the 
surface brightness is low. As was pointed out by Bergvall et 
al. {1999), there is a tendency of a depression relative to an 
exponential light distribution to develop in the centre of the 
galaxies when going to extremely low levels of surface bright- 
ness. Hence, the central disk surface brightness will be inad- 
equate to describe how dim the disk really is. Such cases are 
e.g. ESO 505-04, ESO 546-34 and ESO 602-26. The central 
surface brightness entries have been corrected to face-on pro- 
jection for the purpose of standardization. It should be noted 
that these corrections may be too large for individual objects if 
the disk contains an embedded stellar component with a struc- 
ture different from that of a thin disk, e.g. a bulge. Such effects 
are however negligible for the sample as a whole. 
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Table 1. The integrated properties of our sample of 9 blue LSBGs. The morphological types have been determined by Ronnback 
& Bergvall J1994i . Absolute magnitudes are based on distances derived from the model by Schechter ( 1980 1, assuming a dis- 
tance to the Virgo cluster of 17 Mpc and a Hubble constant of Hq = 75 km s _1 Mpc _1 . Unless stated otherwise, m B represents 
the integrated B-band magnitude inside the fj. B Q - 26.5 mag arcsec~ 2 isophote. Here, represents the "true" central surface 
brightness and /i® the central surface brightness of the disk extrapolated to the centre. Both have been corrected for inclination 
assuming an infinitely thin disk. BVR are in the Cousins system, i in the Thuan-Gunn system and JHK' in the Johnson system. 
All magnitudes have been corrected for galactic extinction (Burstein & Heiles 1982). EW represents the equivalent width of the 
Ha emission line. All photometric data come from Ronnback & Bergvall J 1 9941 and Bergvall et al. ( 19991, except the <B — R> 
entries, which come from the ESO/Uppsala survey (Lauberts & Valentijn 1989 1. 



ESO id. 


Type 


m B 


M B 


(mag. 

arcsec~ 2 ) 


(mag. 
arcsec~ 2 ) 


B—V 


<B-R> 


V-i 


B-J 


J-H 


J-K' 


O/H 


EW 

(A) 


log 

M(HI) 

(M G ) 


074-16 


SBcd 


16.0 


-17.3 


23.9 


23.6 


0.44 


0.42 


0.08 


1.76 


0.5 








9.69 


084-41 


Sd 


17.2 




25.2 


24.7 


0.44 


0.47 


0.08 


1.71 












146-14 


Sdm 


15.4 a 


-16.9 


23.6 


23.8 


0.36 


0.30 


0.02 


1.72 


0.7 


0.7 


7.61 


37 


9.57 


288-48 


Pec. 


16.8 


-14.4 


22.0 


23.3 


0.43 


0.41 


0.02 


1.70 












462-36 


SBd 


16.4 




22.8 


23.1 


0.46 


0.42 


0.10 


1.94 




0.7 








505-04 


Sm 


15.3 


-16.9 


22.7 


21.3 


0.49 


0.42 


-0.01 








7.68 


15 


8.67 


546-09 


Sm 


15.2 


-18.0 


21.9 


22.7 


0.43 


0.43 


0.24 


1.80 






7.95 




9.27 


546-34 


Sm 


15.7 


-15.8 


24.5 


22.8 


0.38 


0.43 


0.01 


1.40 






7.64 


37 


8.81 


602-26 


Pec. 


17.1 




24.7 


23.6 


0.31 


0.39 


0.05 


1.40 


0.7 











a) Magnitude based on the ji B = 25 mag arcsec 2 isophote 



3. Spectral evolutionary modelling 

When analyzing the integrated properties of galaxies, it is im- 
portant to use spectral evolutionary models as closely adapted 
to the known conditions of the target objects as possible. Here, 
the model of Zackrisson et al. (2001 hereafter Z01) is mainly 
used, with the publicly available model PEGASE.2 (Fioc & 
Rocca- Volmerange 1 1 999> employed to a lesser extent to test 
the robustness of all important conclusions. These models are 
adopted because of their capability to predict the spectral en- 
ergy distribution of stellar populations with ages spanning from 
young to very old, as well as the corresponding component of 
nebular emission, which can be crucial for the interpretation of 
optical and near-IR photometry in star-forming galaxies (see 
e.g. Z01, Ostlin et al. l2003l . 

Z01 is based on the method of isomass synthesis and em- 
ploys synthetic stellar atmospheres by Lejeune (1998 1 and 
Clegg & Middlemass (1987), together with stellar evolution- 
ary tracks mainly from the Geneva group. The stellar com- 
ponent includes pre-main sequence evolution and a stochastic 
treatment of horizontal branch morphologies at low metallici- 
ties. Gas continuum and emission lines are predicted using the 
photoionization code Cloudy version 90.05 (Ferland 1996 1. For 
each time step, the spectral energy distribution of the integrated 
stellar population is used as input to Cloudy, which makes the 
computation time-consuming, but provides realistic temporal 
evolution of the nebular component. 

PEGASE.2 is based on the method of isochrone synthe- 
sis, employs the same synthetic stellar atmospheres as Z01, 
but stellar evolutionary tracks mainly from the Padova group. 
Using a simple prescription for the nebular component, it pre- 



dicts both continuum and line emission of the ionized gas, al- 
though with non-evolving emission line ratios. 

3.1. Extinction corrections 

We apply corrections for intrinsic extinction by assuming an 
extinction curve characteristic of the Milky Way (Miller & 
Mathews 1972 1. Since the extinction, based on the measured 
Balmer decrement (Ronnback 1993 1, appears to be very low in 
these objects, we consider c = log w (Fu£fl/ F^b) = 0.15 + 0.05 
(corresponding to E(B -V) - 0.03-0.07 mag), where F^p rep- 
resents the observed HB emission line flux and Fjjp,o the HB 
emission line flux in the absence of extinction. 

3.2. Chemical abundances 

Because of the very low metallicities of blue LSBGs (5-20% 
solar) and the small metallicity dispersion among their dif- 
ferent HH-regions (Ronnback 1993), constant metallicities of 
Z = 0.001 and Z = 0.004 are adopted for the grid of evolu- 
tionary sequences used. A more realistic approach would be to 
use a self-consistent evolution of the metallicity, and allow the 
rates of infall of primordial gas or blow-out of metal-enriched 
material to be free parameters. The measured metallicities for 
each galaxy could then be used to determine the quality of a 
given fit. This would not give rise to a constant metallicity, but 
rather a metallicity distribution with a luminosity-weighted av- 
erage corresponding to the measured value. Such an approach 
would however increase the required grid of model evolution- 
ary sequences by a large factor and is not essential here. 
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3.3. The initial mass function 

The initial mass function (IMF) of the stellar component is as- 
sumed to follow a single-valued power law 



<P(m, a) — — am 
dm 



(1) 



inside the mass interval [M\ ov/ , M up ], where a slope of a = 
2.35 corresponds to the commonly adopted Salpeter IMF. The 
Salpeter slope is often assumed to provide an adequate descrip- 
tion of stars more massive than « 0.5 M for both the Galaxy 
and the Large Magellanic Cloud, but has been shown to signif- 
icantly overpredict the number of stars less massive than this 
(e.g. Kroupa 2001 1. Although the necessary modification of the 
IMF at the lowest masses has a pronounced effect on the mass 
budget of normal stellar populations, the affected low-mass 
stars contribute very little to the integrated light. A common 
strategy is therefore to use a single-valued power law when 
calculating the spectral energy distribution of synthetic stellar 
populations, and to rescale the mass-dependent predicted quan- 
tities (e.g. stellar M/L ratios) to more realistic IMFs afterwards. 
This is the procedure adopted in the following. 

It should however be noted that much probably remains to 
be learned about the IMF and its cosmological evolution. It 
has been argued that the IMF could have been more top-heavy 
in the early universe (e.g. Larson [T998i . with a slope flatter 
than the Salpeter value by Aa « 0.4. Conversely, it has also 
been pointed out that the local high-mass IMF may actually be 
steeper than the Salpeter value by Aa a 0.4, due to unresolved 
binary systems (e.g. Scalo ll998l Kroupa l200H . To test the ef- 
fects of IMF variations within this range of plausible slopes, 
spectral evolutionary scenarios based on single-valued IMFs 
with a = 1.85-2.85 are therefore used in this study. 

Although the upper mass limit M up is often assumed to lie 
around 120 M , the precise value and its cosmic variance re- 
main very uncertain due to small number statistics. For this 
reason, spectral evolutionary sequences with M up = 40 M are 
also tested. 



3.4. Star formation history 

Three different kinds of star formation histories (SFHs) are 
considered: a star formation rate (SFR) which stays constant 
during a given time and then drops to zero; an exponentially de- 
creasing SFR (<x exp(-f/r), with e-folding decay rate t > 0); 
and an exponentially increasing SFR (same, but with t < 0). 
The first is suitable for abrupt starbursts (with duration of bursts 
typically 10-100 Myr), whereas the second is more appropriate 
for late-type galaxies with signs of star formation over signif- 
icant fractions of the age of the universe. The third has been 
inspired by the models of Boissier et al. (I2003i in which low 
surface brightness galaxies differ from high surface brightness 
ones in the larger angular momenta of their dark matter halos. 
Their models predicts that the SFRs of the galaxies with low- 
est mass and highest spin should increase by several orders of 
magnitude over a Hubble time. Bell et al. (2000 1 also find evi- 
dence that an increasing SFR is required to explain the colours 
of certain LSBGs. Although the SFHs predicted in Boissier et 



Table 2. The grid of Z01 evolutionary sequences. All se- 
quences assume Z gas = Z stms , a gas covering factor of 1.0, a 
hydrogen number density of 100 crrr 3 and the IMF slope, a, 
to be valid throughout the stellar mass range 0.08-120 M Q . 
The grid consists of all possible combinations of the parameter 
values listed below. In addition to these, 108 sequences with 
M up = 40 M Q were generated for test purposes with variations 
in IMF slope, metallicity, star formation history and redshift, 
but with fixed M tot = 10 9 M Q and filling factor 0.01. 



M tol (M s ) 
IMF, a 
SFH 



10 9 , 10 10 
1.85, 2.35, 2.85 

clO 8 , e5 x 10 s , el x 10 9 , e3 x 10', e6 x 10' 
el5 x 10", cl5 x 10', e-15 x 10 9 , e-2 x 10 9 
Metallicity, Z 0.001, 0.004 
Filling factor 0.01,0.1 
Redshift, z 0.00, 0.01 
IMF: dN/dm cc mT" 
SFH=Star formation history. 

c=Constant star formation rate during the subsequent number of 
years. 

e=Exponentially declining/increasing star formation rate with t equal 
to the subsequent number of years. 



Table 3. The grid of PEGASE.2 evolutionary sequences. All 
sequences assume the IMF slope, a, to be valid throughout the 
stellar mass range 0.08-120 M . The grid consists of all pos- 
sible combinations of the parameter values listed below. The 
abbreviations used are described in Table |2| 

IMF, a 1.85, 2.35, 2.85 

SFH clO 8 , e5 x 10 8 , el x 10 9 , e3 x 10 9 , e6 x 10 9 

el5 x 10 9 , cl5 x 10 9 , e-15 x 10 9 , e-2 x 10 9 
Metallicity, Z 0.001, 0.004 



al. J2003t cannot be perfectly reproduced using an exponen- 
tially increasing SFR, we find that a t = -2 Gyr scenario 
increases the SFR by roughly the right amount over the age 
of the universe to accord with their steepest SFH predictions. 
Scenarios with a constant SFR or an exponentially increasing 
one with r = — 15 Gyr are considered adequate to reproduce 
their less dramatic SFH slopes, except during the earliest phase 
of evolution. 



3.5. Other model parameters 

In Z01, the properties of the nebular component is regulated 
by several parameters not considered in PEGASE.2, some of 
which can have a significant effect on predicted broadband 
magnitudes (see Z01 for a demonstration): the metallicity of 
the gas, the covering factor, the filling factor, the gas density 
and the total mass of the stellar component. Here, the metal- 
licity of the g assumed to be the same as that of 
the stellar population, Z stars . Since stars pollute the interstellar 
medium with metal-enriched material, models of chemical evo- 
lution generally predict the metallicity of the stars to be lower 
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than that of the surrounding gas. The approximation used here 
is nonetheless adequate as long as the metallicity is low. 

The covering factor regulates the fraction of the stellar com- 
ponent covered by gas (in terms of surface area). Since blue 
LSBGs are gas-rich and no studies to date have discovered 
Lyman-continuum leakage from any LSBG, it is likely that 
the covering factor is very high. Here, we therefore consider 
a cover factor of unity only. 

The filling factor, defined as the fraction of the total vol- 
ume of the Stromgren sphere occupied by gas, is constrained by 
observations to lie in the range 0.01-0.1 (Martin H 9971 Miller 
1996 1, which is the interval adopted here. 

The hydrogen number density of the gas condensations in- 
side the Stromgren sphere, has been taken to be 100 cirT 3 , 
based on Copetti et al. d2"0"0lfl . 

In Z01, the nebular component is computed under the as- 
sumption that the stars of a galaxy are concentrated to the cen- 
tre of a single spherical nebula, rather than distributed among 
a large number of smaller Stromgren spheres. This approxi- 
mation gives reasonable results for the hydrogen Balmer lines 
used in this study, but may not be as good for emission lines 
more sensitive to the detailed temperature and ionization struc- 
ture of the cloud. In the framework of this approximation, the 
number of Lyman continuum photons striking the illuminated 
face of the cloud at a given population age, t, is regulated 
by SFR(f) and the total gas mass, M tot , which may eventu- 
ally be converted into stars during the course of a given SFH. 
Since previous investigations have found the baryonic mass of 
blue LSBGs to he in the interval 10 9 -10 10 M Q (De Blok & 
McGaugh ll997> . this is the range adopted here. In the case of 
star formation rates increasing with time, M tot is defined as the 
total gas mass that would have been consumed after 15 Gyr. It 
should be noted that Z01 does not assume that gas expelled 
from the stellar population can be recycled into subsequent 
generations of stars. 

Finally, the conversion from synthetic population spectra to 
fluxes in broadband filters depends on the redshift of the target 
galaxies. Due to the galaxy diameter limit of the ESO-Uppsala 
catalogue, none of our objects are expected to lie significantly 
further away than z ~ 0.01. For this reason, only redshifts in 
the range z = 0-0.01 are considered. 

3.6. Model grids 

Using Z01, a grid of 540 evolutionary sequences has been gen- 
erated with parameter values summarized in Tabled constitut- 
ing the most extensive set of spectral synthesis predictions used 
for the study of LSBGs so far. In addition, a smaller set of 54 
evolutionary sequences, with parameter values summarized in 
Table |3] was computed using PEGASE.2 to be used for con- 
sistency testing. Each sequence runs from ages of 1 Myr to 15 
Gyr in small time steps. 

3. 7. Ha equivalent widths 

In addition to broadband fluxes, EW(Ho-) are also used as a 
constraint when determining the properties of the stellar pop- 



ulations of the galaxies in our sample. Although the required 
Ha narrowband data are available only for three of our nine 
objects (with EW(Ha) = 37, 15 and 37 A respectively), the 
global EW(Hof) of blue LSBGs appear to typically lie inside 
a limited interval. Narrowband images of two additional blue 
LSBGs included in the original sample but lacking photometry 
(ESO 359-3 1 and ESO 405-06), reveal EW(Ha) of 19 and 25 A 
respectively. The equivalent widths derived from VLT/FORS2 
long-slit spectra along the major axis of four highly inclined 
LSBGs (ESO 031-13, ESO 462-32, ESO 532-32, ESO 548- 
04) with colours very similar to those in the present sam- 
ple (Zackrisson et al., in preparation) are EW(Ha) = 30, 39, 
43 and 18 A. Due to the largely transparent disks of LSBGs 
(e.g. Matthews & Wood 2001 1, these long-slit observations 
should provide reasonable estimates of the global EW(Har). 
With o-(EW(Ho')) = 20%, the EW(Ha) measurements for aU 
9 galaxies approximately lie inside the interval 10-60 A at the 
2cr level, which is the adopted range of EW(Ha) considered 
acceptable for successful model fits in this study. 

3.8. Homogeneity of the stellar population 

The small colour gradients in our objects indicate that their 
stellar populations must be fairly homogeneous. The small dis- 
persion in EW(Ha) among blue LSBGs furthermore shows 
that their star formation history cannot have proceeded very 
stochastically, i.e. through violent starbursts, unless the IMF is 
peculiar. 

However, in a study of blue compact galaxies, Bergvall 
and Ostlin ( 2002 1 detected abnormally red colours, typical of a 
metal-rich population, in the faint halos of these objects (which 
are generally considered to be very metal-poor). One remote 
possibility is that blue LSBGs may have formed recently as a 
consequence of a merger between a relatively metal-rich early- 
type galaxy and an intergalactic gas cloud. In such a scenario, 
the stars formed in the merger event could very well be outshin- 
ing the older population, giving the appearence of a genuinely 
young galaxy. The old, underlying host would then only be ob- 
servable in the very faint outskirts of the blue LSBG, analo- 
gous to the case of blue compact galaxies. This scenario may 
however be tested by studying the halo colours of nearly edge- 
on galaxies, of which ESO 146-14 from our sample is one ex- 
ample. Within the error bars, we do not detect any significant 
gradients in the optical/near-IR colour profiles in the direction 
perpendicular to the major axis of the disk. For this reason, 
mixed populations more extreme than those already produced 
by the adopted range of star formation histories, e.g. by adding 
a single-age young population to a single-age old one, are not 
considered in the present study. 

3.9. Model fitting 

For each time step along each of the 540 Z01 evolutionary se- 
quences defined in Table [2] and the 54 PEGASE.2 sequences 
of Table a least-squares fit is performed for each object be- 
tween the N observational and model passband fluxes, / b s ,, and 
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Fig. 1. The temporal evolution of B — i vs. V - J predicted by Z01 (solid) and PEGASE.2 (dashed). Both models assume a 
metallicity of Z = 0.001, a Salpeter IMF throughout the mass range 0.08-120 M Q , a redshift of z = and an exponential star 
formation history with a) t = 1 Gyr, b) t = 6 Gyr and c) t — -2 Gyr. The additional parameter values for Z01 are a filling factor 
0. 1 and a total gas mass of 10 10 M available for star formation. Markers indicate ages of 1 Myr (circle), 100 Myr (triangle), 1 Gyr 
(square), 5 Gyr (pentagram, omitted from rightmost plot to avoid cluttering) and 15 Gyr (asterisk) along the model sequences. 
The colours of the eight objects (corrected for intrinsic extinction using c = 0.15), for which data exist in all four filters required, 
are indicated by crosses representing lcr error bars. 



/mod,;, weighted by their photometric uncertainties, cr ; . An aver- 
age model-observation residual, A mo( j, is then calculated using: 

Amod " V^VT' (2) 

where w>; = crj 2 and 

. /./obs,/ — ./mod,/ \ 



To estimate confidence levels on the stellar population pa- 
rameters derived from this fitting procedure, numerically esti- 
mated lcr (68%) and 2cr (95%) levels on the model-observation 
residuals are imposed. 

During the fitting procedure, the observed photometric 
fluxes in all filters except R are used. The reason for the omis- 
sion of /?-band data is the uncertain calibration of low surface 
brightness objects in the ESO/Uppsala catalogue. When com- 
paring the B - R photometry from the ESO/Uppsala catalogue 
with the results from CCD photometry of bulge-dominated 
LSBGs carried out by Beijersbergen et al. (1999 1, a rather 
large scatter in the residuals becomes evident, corresponding 
to <tb-r ~ 0.3 mag. For our galaxies which have fainter central 
regions, the calibration problems may be even more severe. 

Although not explicitly used in the fitting procedure, Ha 
data are used to reject all fits with emission line equivalent 
widths outside the typically observed range of 10 A < EW(Ha) 
< 60 A, as discussed in section 3.7. To allow for the uncertainty 
in the extinction correction, values of c = 0.10, c = 0.15 and 
c = 0.20 are considered for each galaxy. 

4. Results 

In Fig.^i and[Q>, we present the behaviour predicted by Z01 
and PEGASE.2 in a diagram of B - i vs. V - J for two of the 
evolutionary scenarios most successful in reproducing the ob- 
served colours of the galaxies under the assumption of a stan- 



dard Salpeter IMF. The blue LSBGs in our sample are all lo- 
cated inside a small region of this two-colour diagram, indicat- 
ing that they all share similar properties. Although slight differ- 
ences between the model predictions of Z01 and PEGASE.2 do 
exist in these filters, both models are able to roughly reproduce 
the observed colours of our objects. 

The large number of spectral evolutionary sequences con- 
sidered in our analysis allows several hundreds to thousands 
of fits of comparable quality. Among these, a wide variety of 
model parameter values can be found, indicating that the ob- 
servations considered here can be equally well reproduced by 
many fundamentally different star formation scenarios. This 
degeneracy becomes obvious in Fig. and \T]p, where the 
colours of our target objects may be reproduced with both 
star formation histories depicted (exponentially decaying SFRs 
with t — 1 and 6 Gyr respectively), although at different ages 
along these evolutionary sequences. In the r = 1 Gyr case 
(Fig.^), all objects appear to have ages somewhere between 1 
and 5 Gyr, whereas the t = 6 Gyr case (Fig.^) implies that all 
objects should have ages in excess of 5 Gyr. This degeneracy 
between age and star formation history is an old problem (see 
e.g. McGaugh & de Blok lT997l Gil de Paz & Madore l20T)2l . 
which nonetheless too often tends to be overlooked. The inclu- 
sion of near-IR broadband data, which should be a better tracer 
of old stars, does not appear to improve the situation substan- 
tially. A similar situation also occurs for parameters other than 
the star formation history, and as it turns out, neither the slope 
nor the upper mass limit of the IMF can be constrained by the 
available observations. 



4.1. Star formation history 

Although the available photometry does not appear to impose 
any strong constraints on the star formation history of the blue 
LSBGs, their low EW(Ha) do. In Fig. [2 we compare the tem- 
poral evolution of EW(Ha) predicted by Z01 and PEGASE.2 
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Fig. 2. The temporal evolution of EW(Ho') predicted by Z01 (solid lines) and PEGASE.2 (dashed) for various SFHs compared 
to the range observed among the bluest LSBGs (gray area). Both models assume Z = 0.001. Z01 furthermore assumes a filling 
factor of 0.1 and M tot = 10 10 M@. The three panels represent different assumptions about the IMF. a) A Salpeter IMF with 
M up = 120 M Q . b) An IMF slope of a = 2.85 with M up = 120 M Q . c) A Salpeter IMF with an upper mass limit of 40 M Q . 
From top to bottom in the right half of each panel, the different lines indicate star formation histories with r = — 2 Gyr, -15 Gyr, 
constant SFR, r = 15 Gyr, 6 Gyr, 3 Gyr, 1 Gyr and 500 Myr. The line dropping rapidly at an age of 100 Myr corresponds to the 
short burst scenario. The dash-dotted line running up to 1 Gyr in a) represents the predictions of Starburst99 for a constant SFR. 



for different SFHs, with the range of EW(Hor) observed among 
blue LSBGs. 

In the short burst scenario (constant SFR for 100 Myr, fol- 
lowed by no star formation thereafter), much too high values of 
EW(Ha) are predicted during the early phase of active star for- 
mation. Once star formation has ended, EW(Ho') on the other 
hand rapidly drops to values much too low, allowing a time 
span of only a 20 Myr during which this scenario predicts 
equivalent widths in the observed range. However, even with 
the wide range of IMFs and metallicities considered in this pa- 
per, the bluest LSBGs are not blue enough (in e.g. B - J) to 
be consistent with the colours predicted by the models for this 
intermediate phase. Hence, the short burst scenario (as defined 
here, with star formation ending very abruptly) can be ruled 
out. 

In scenarios including constant or increasing star forma- 
tion rates over cosmological time scales, such as those moti- 
vated by predictions for the lowest-mass low surface brightness 
galaxies in the model of Boissier et al. (|2003 1, the EW(Ha) 
never have time to drop sufficiently to reach the observed range 
within the age of the universe (14. 1^9 Gyr; Tegmark et al. 
2004 1 in the case of an IMF with a Salpeter or flatter (e.g. 
a = 1.85) slope. For the most rapidly increasing star formation 
rate (r = -2 Gyr) the same IMFs furthermore produce colours 
which are much too blue for our objects, as demonstrated for a 
Salpeter IMF in Fig.[0;. Both these conclusions are confirmed 
by PEGASE.2, indicating that, under the assumption of a stan- 
dard IMF, these star formation scenarios cannot be valid for 
the low-metallicity, low surface brightness galaxies considered 
here. The idea of a highly variable past SFH, which was ad- 
vanced by Boissier et al. (2003 1 to explain the presence of red 
LSBGs, does not appear attractive for our objects, as discussed 
in Section 3.8. In Fig. [2^, the EW(Hff) evolution during the 
first Gyr predicted by Starburst 99 (Leitherer et al. 19993 ver ~ 
sion 4.0), for a constant SFR scenario, is included to demon- 
strate that the problem is not due to unrealistically high Lyman 
continuum fluxes in Z01 and PEGASE.2, which employ stel- 



lar atmosphere models for high-mass stars considerably less 
sophisticated than those of Starburst 99. If anything, Z01 and 
PEGASE.2 appear to underpredict the equivalent widths. 

Although the EW(Ha) predicted by Z01 and PEGASE.2 
neglect the effect of underlying absorption by the stellar popu- 
lation, this may have increased the predictions by no more than 
~ 5 A (Ronnback & Bergvall 1995 1 compared to the obser- 
vations for these objects. Selective extinction due to emission 
lines originating in regions more dust-rich than those respon- 
sible for the stellar continuum (Calzetti et al. 1994 1 may fur- 
thermore lead to corrections of only < 25% to the observed 
EW(Ha) for the low levels of extinction (Ronnback 1993 1 ob- 
served in our targets. 

The very high EW(Ha) predicted by scenarios including 
constant or increasing star formation rates may nonetheless be 
remedied if the slope of the IMF is significantly steeper (e.g. 
a = 2.85), or the upper mass limit substantially lower (e.g. 
M up - 40 Mo) than typically assumed, as shown in Fig. I2J3 
and|2;, respectively. Lowering M up does not necessarily contra- 
dict the observed (N/O) and (O/H) abundances of blue LSBGs 
(Ronnback 1993 1, since this alteration has only a modest im- 
pact on the abundances of these elements (Olofsson 1993. A 
bottom-heavy (high a) IMF could help to explain the high M/L 
of LSBG disks claimed by Fuchs ( 2002 1, as demonstrated by 
Lee et al. (2004), but may produce too low (O/H) abundances. 
A high a just for high-mass (> 1M Q ) stars would on the other 
hand be consistent with recent models for the field-star IMF 
(Kroupa & Weidner 2003). A more detailed investigation of 
what constraints the EW(Ho') and chemical abundances can im- 
pose on the SFHs and IMFs of LSBGs (and not just the bluest 
objects) would be very valuable. 

An alternative way to reconcile the median observed 
EW(Ha) with constant or increasing SFRs could be to reduce 
the hydrogen-ionizing flux shortward of the Lyman limit by ad- 
vocating dust extinction. In order for this to work, the required 
Lyman continuum depletion factor would however have to be 
higher than 2.5. It is not obvious that the amount and spatial 
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distribution of dust required to achieve this is at all realistic 
for blue LSBGs, whose star-forming regions show only mod- 
est signs of reddening. 

4.2. Ages 

Here, the absolute age (or age, for short) of a galaxy is defined 
as the time passed from the first onset of significant star forma- 
tion processes, i.e. the start of the assumed SFH. The average 
age of a galaxy, which will be discussed in relation to previous 
investigations on LSBGs, is instead defined as the average age 
of all stars ever formed in the system. Since stellar remnants 
are also included in this inventory, the average age is a function 
of SFH and absolute age only. 

The range of absolute ages among the acceptable fits gen- 
erated by Z01 and PEGASE.2 for our blue LSBGs are pre- 
sented in Table 0] These intervals all assume a Salpeter IMF 
with M up = 120 M Q . 

Both models agree that the observations for most of the ob- 
jects in our sample can be equally well reproduced by ages 
anywhere in the range from 1 or 2 Gyr, up to 15 Gyr, which 
is the highest age considered here. The ages are therefore not 
well constrained by the available observations, despite the use 
of near-IR data. 

What about the average ages? In Bell et al. (2000), colour- 
colour diagrams are presented with model grids containing 
variations in metallicity and star formation history. For an as- 
sumed absolute age of the galaxies (12 Gyr, in their case), they 
show that the average stellar age can in fact be constrained to 
impressive levels using broadband photometry in four filters 
only. As we demonstrate in Fig. [2] problems do however arise 
when the absolute age is no longer assumed to be known. In this 
case, the grid shifts as a function of absolute age, so that e.g. 
the colours of the group of blue LSBGs around B - V » 0.4 
mag, V — J ~ 1.2 mag in Fig. [5] may be equally well repre- 
sented by an absolute age of 3.0 Gyr, an average age of 2.2 
Gyr and a star formation history with r = 1 .0 Gyr (left panel) 
as by an absolute age of 13.0 Gyr, an average age of 7.4 Gyr 
and t = 15.0 Gyr (right panel). Hence, neither the absolute nor 
average ages are particularly well constrained. For the LSBGs 
studied by Bell et al. ( 2000), the uncertainty is likely to be even 
more severe due to the unknown metallicities of their objects. 

From the number of missing entries in Tabled it is evident 
that both Z01 and PEGASE.2 have some trouble explaining 
the observed properties of certain galaxies even at the 2cr level, 
although the problematic objects are different for the two mod- 
els. Both Z01 and PEGASE.2 identify about three objects each 
which are better fitted by young ages (< 3 Gyr), but these are 
also the most difficult galaxies to fit simultaneously with both 
codes. Therefore, the two models cannot be said to convinc- 
ingly confirm each other's detections of young galaxies. 

4.3. Stellar mass-to-light ratios and population masses 

The mass-to-light ratios (M/L F ) used here are defined as 
W I Z^' wnere F represents a particular filter and M the mass 
locked up in stars or stellar remnants. From the V-band mass- 



to-light ratios (M/Ly,) and the V-band luminosities (Ly) pre- 
dicted by the model fits, stellar population masses (M*) are 
computed. Table |4] lists the range of M/Ly and M* allowed 
by the acceptable fits generated by Z01 and PEGASE.2 for our 
target objects. At the 2<x level, M/Ly = 0.5-3.5 and M* * 10 7 - 
10 10 M Q . These constraints all assume a Salpeter IMF with 
M up = 120 M . 

The derived M/L ratios are of course subject to substan- 
tial uncertainties due to the poorly constrained scaling of 
the IMF. A single-valued Salpeter power-law is generally be- 
lieved to produce too many low-mass stars (e.g. Kroupa 12001 1 
Chabrier 2003 1, and scaling relations of the type (M/L) w 
X (M/L)sai P eter have been suggested. Both direct integration of 
the observed mass function (Kroupa 2001 ) and dynamical argu- 
ments (McGaugh 2004 1 favour a value of X » 0.6 (factor relat- 
ing a Salpeter IMF with a lower mass limit of M\ ovl = 0.08 M , 
as adopted here, to the Kroupa 2001 IMF), and it may therefore 
be suitable to multiply the M/L ratios of Table [5]by a factor of 
this order. 

In many contexts, e.g. in rotation-curve decompositions and 
in the study of the Tully-Fisher relation, it is advantageous to 
use the filter passband for which the stellar M/L ratios are the 
least vulnerable to uncertainties in age, dust, metallicity and 
star formation history. In de Jong ( 1996) and Bell & de Jong 
(2001} it was argued that the A'-band should be superior in 
this respect. In Table [5] we present the range of M/L in fil- 
ters BVRHJHK allowed by Z01 and PEGASE.2 for our tar- 
get objects. Due to differences in the stellar evolutionary tracks 
adopted by the two codes (discrepancies between Geneva and 
Padova tracks, compounded by the non-standard treatment of 
horizontal branch morphologies used in Z01), Z01 is somewhat 
underluminous at ages above 1 Gyr and systematically predicts 
higher M/L ratios. 

Although we do find that the allowed range of M/L typi- 
cally becomes smaller with increasing central filter wavelength 
for each spectral evolutionary model separately, this advantage 
is counterbalanced by the systematic differences in the M/L- 
predictions between Z01 and PEGASE.2. The ratio of highest 
to lowest allowed M/L among both models instead has a min- 
imum around the i or /-band. This indicates that while stel- 
lar populations theoretically are expected to display a narrow 
range of near-IR M/L ratios, this effect is partly cancelled by 
model uncertainties associated with the treatment of late stages 
of stellar evolution, which have a significant impact on the 
spectral energy distribution in the near-IR for old stellar popu- 
lations (see e.g. Portinari 2004). A similar conclusion about the 
superiority of the / band in stellar-population M/L studies has 
previously been reached by Worthey ( 1 994 ), by considering the 
age and metallicity sensitivity of M/L in different passbands. 

4.4. The effects of a peculiar IMF 

What happens to the derived quantities if the assumption of a 
Salpeter IMF with M up = 120 M is relaxed? In Fig.0 all the 
combinations of age and M/Ly allowed for the blue LSBGs at 
the lcr level by Z01 as well as PEGASE.2 are plotted, under 
various assumptions about the IMF. 
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Table 4. The range in age, M/Ly and M+ allowed by the acceptable fits to the observations at the lcr and 2<x levels using the Z01 
grid of evolutionary sequences defined in Tableland the PEGASE.2 grid of evolutionary sequences defined in Table|3] Only the 
results from sequences assuming a Salpeter IMF with M up = 120 M G have been included. Missing entries indicate the lack of 
acceptable fits at a given confidence level or, in the case of M*, the lack of redshift data for a particular object. To obtain values 
of M/Ly and consistent with the observed IMF of Kroupa (2001 1, these values should be multiplied by a factor of » 0.6, as 
discussed in section 4.3. 
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a) 3a confidence levels used, due to lack of acceptable fits within the 2a limit 



Table 5. The allowed range of stellar M/L ratios among the blue LSBGs in our sample for which the Z01 and PEGASE.2 models 
are able to provide adequate fits at levels of lcr and 2<x, respectively. Only the results from model sequences assuming a Salpeter 
IMF with M up = 120 M Q have been used. To obtain values of M/L consistent with the IMF of Kroupa ( 2001 ), these values should 
be multiplied by a factor of * 0.6, as discussed in section 4.3. 
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la 
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When keeping M ap = 12OM , but allowing the slope of 
the IMF to deviate from the Salpeter value, the derived span 
of M/Ly is increased upwards, allowing M/Lv as high as 6.1 
(Z01) and 3.6 (PEGASE.2) at the lcr level. The change is al- 
most entirely due to the inclusion of bottom-heavy (a = 2.85) 
IMFs. In this scenario, high M/Ly arise naturally because of 
the increased number of low-mass, low-luminosity stars. Since 
SFHs with constant or increasing SFRs produce acceptable 
EW(Ha) in this case (see Fig. 13, low M/Ly at high ages are 
however also allowed. 



Assuming an upper mass limit below 120 M alters the 
shape of the allowed region in the plot of age versus M/Ly, 
but has a modest impact on the extremum values. Analogous 
to the case with a = 2.85, low M/Ly are produced also at high 
ages. Since the number of low-mass stars is almost unaffected 
in this case, no high M/Ly are however predicted. 



The most notable effect of simultaneously setting M up = 
40 M and allowing a = 1.85-2.85 is that M/Ly as low as 
M/Ly = 0.44 are generated. 



4.5. Comparison to other investigations 

Despite the use of both optical and near-IR broadband photom- 
etry, the age constraints derived here (2-15 Gyr) are essentially 
equivalent to those derived mainly from optical photometry by 
Ronnback ( 1993 1 and Bergvall & Ronnback ( 1994 1, indicating 
that the addition of near-IR photometry is inadequate to impose 
useful constraints on the ages without prior knowledge of the 
star formation history. The mass estimates for the stellar com- 
ponent presented here are also of the same order of magnitude 
as those derived in Bergvall & Ronnback ( 1994 1 and Ronnback 
( 1993 ), although estimates vary significantly for individual ob- 
jects due to some of the very extreme IMFs allowed in these 
older papers. 

The majority of LSBGs adhere to the Tully-Fisher relation, 
indicating that their total M/L ratios are higher than in normal 
disk galaxies and that these objects are more dark matter dom- 
inated (e.g. McGaugh & de Blok 1998). This feature has made 
LSBGs popular targets for measuring the dark halo density pro- 
file. In the procedure of rotation curve decomposition common 
when trying to distinguish between different density profiles, 
some M/L ratio believed to be representative of the stellar com- 
ponent is typically assumed. Often, this ratio is adopted from 
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Fig. 3. Diagrams of B - V vs. V - J, showing the predicted positions of stellar populations at absolute ages of 3 Gyr (left panel) 
and 13 Gyr (right panel). In each panel, lines corresponding to metallicities of Z = 0.001 (thick solid) and Z = 0.004 (thin solid) 
connect the colours predicted for different star formation histories, parametrized by r in SFR(f) oc exp (— f /t). Each combination of 
absolute age and t corresponds to the listed average stellar population ages, <Age>. Dashed lines connect predicted populations 
with the same SFH but different Z. The colours (corrected for intrinsic extinction using c = 0.15) of the eight objects for which 
data exists in all three filters are indicated by crosses representing lcr error bars. A Salpeter IMF with M up = 120 M is assumed. 



colour-M/L relations (e.g. Bell & de Jong 2001 1, which are 
based on spectral evolutionary scenarios. The M/L ratios de- 
rived this way do however depend on the assumptions regard- 
ing the properties of the stellar populations of LSBGs, such as 
the age, metallicity and star formation history used when con- 
structing the relation. 

When the colour-M/L relations from Bell & de Jong ( 200 1 ) 
are applied to our targets (B — V = 0.24-0.46 mag after cor- 
rection for intrinsic extinction c = 0.1-0.2), stellar M/L of 
M/Ly = 0.5-1.6 and M/L K = 0.5-0.8 (from their Table 3, 
after correcting by a factor of 1 .57 for differences in the scaling 
of the Salpeter IMF) are predicted, which is consistent with the 
best-fit M/L v = 0.5-3.1 and M/L K = 0.4-2.0 derived here. 
Our range of allowed M/L is more generous, but this is ex- 
pected, since these Bell & de Jong predictions only take into 
account the scatter among average colour-M/L relations de- 
rived from different evolutionary scenarios, whereas our pre- 
dictions include the total scatter due to the fitting procedure as 
well as systematic differences between two spectral evolution- 
ary models. There is furthermore no indication that the M/Lk 
values derived by us are in conflict with the maximum disk val- 
ues (Fig. 6 in Bell & de Jong 2001 ) inferred from observations 
of LSBGs at the bluest colours. 

5. Discussion 

Although the available data appear insufficient to claim a young 
age for the blue LSBG population, ages of 2-5 Gyr are cer- 
tainly allowed for all the galaxies in our sample and possibly 
favoured for a few objects. In the framework of an Q.m = 0.3, 
Qa = 0.7 (ACDM) universe, an age this low would correspond 
to a formation redshift of z < 0.5. According to the hydro- 



dynamic galaxy formation models of Nagamine et al. l|2001 1, 
only a very small fraction of galaxies with present mass in stars 
lower than 10 10 M Q (such as the blue LSBGs in our study) 
should form at such low redshifts in this cosmology. 

To estimate the severity of this potential problem, the frac- 
tion of blue LSBGs in the local universe should be compared 
to the fraction of objects with similar masses formed at low 
redshift in the hierarchical models. 

From the ESO-Uppsala catalogue (Lauberts & Valentijn 
1989), one may estimate that the blue LSBGs obtained with 
our selection criteria constitute « 18% by number per unit vol- 
ume of those with Hubble types from Sa to Im. From a sample 
of 1550 galaxies with M B < -12.5 in nearby galaxy clusters, 
Ferguson & Sandage ( 1991 1 estimate that the galaxies of these 
Hubble types constitute w 20% of the total galaxy population 
(most of the rest being dwarf ellipticals). This would indicate 
that the local fraction of blue LSBGs is at least as 3.2% by num- 
ber. Due to the underrepresentation of S a to Im types in clusters 
compared to lower-density environments, where many LSBGs 
may be located (Rosenbaum & Bomans 2004 1, this fraction 
should be conservative. 

Assuming B-V ~ 0.7 mag (Karick et al. 2003 1 and a stellar 
population M/L v sa 1.5 (Mateo et al. 119981 for dE/dSph dwarf 
galaxies, the Ferguson & Sandage luminosity limit corresponds 
to a stellar mass of « 2.4x 10 7 M . For the smaller fraction of dl 
galaxies, an HI mass of » 2 x 1 7 M Q (Hoffman et al. 1 996 ) may 
be inferred at this luminosity. Given typical gas mass fractions 
for dl galaxies (0.4-0.8; Pilyugin & Ferrini 2000 1, this converts 
into stellar population masses not very different from those of 
dE/dSphs. By taking the risk of overinterpreting small number 
statistics, one may then from Fig. 4 in Nagamine et al. (2001 1 
estimate that only « 0.2% of the objects with stellar masses 
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Fig. 4. The allowed ages and M/Ly of acceptable fits (markers) for both Z01 and PEGASE.2 at the lcr level, under various 
assumption about the IMF: a) M up = 120 M Q , a = 2.35 for Z01, b) M up = 120 M Q , a = 1.85-2.85 for Z01, c) M up = 40 M , 
a = 2.35 for Z01, d) M up = 40 M @ , a = 1.85-2.85 for Z01, e) M up = 120 M , a = 2.35 for PEGASE.2, f) M up = 120 M s , 
a = 1.85-2.85 for PEGASE.2. 



2.4 x 10 7 < M/M Q < 10 10 should form at z < 0.5 (correspond- 
ing to an age < 5 Gyr) and « 0.8% at z < 1.0 (age < 7.5 Gyr). 
The typical redshift of formation for objects of these masses is 
instead predicted to be z ~ 3-5.5 (corresponding to ages of 1 1- 
12 Gyr). Actually, the formation redshifts used in Nagamine et 
al. (|2001 1 correspond to the mass-weighted averages of the for- 
mation epochs of the stellar particles in their simulation, mak- 
ing this comparison to our absolute age estimates (which are 
by definition systematically higher) very conservative. 

We therefore conclude that if the blue LSBGs in our sam- 
ple really have ages below 5 Gyr, there would be an order of 
magnitude (3.2% compared to 0.2%) too many such objects in 
the local universe to be compatible with the Nagamine et al. 
(2001 1 predictions. If our targets are younger than 7.5 Gyr, the 
discrepancy is instead a factor of 4 (3.2% compared to 0.8%). 
Hence, the ages of the bluest LSBGs may potentially provide 
an important test of galaxy formation scenarios. Since LSBGs 
typically are non-merging disk galaxies with well-behaved ro- 
tation curves, they are likely to be superior in this respect com- 
pared to other young galaxy candidates, like blue compact and 
Ffll-galaxies, for which estimates of both stellar and total mass 
may be more problematic. Developing methods to constrain the 
age of the blue LSBG population with higher precision would 
therefore be very valuable. 

What observations would be required to better constrain the 
ages? The strength of the 4000 A break, D4000, is often used 
as an age indicator for stellar populations, but actually depends 
very weakly on age at low metallicities (e.g. Worthev 1 19941 



Maraston 2004 1. Hydrogen Balmer lines in absorption appear 
potentially useful (e.g. Worthey 1994, Vazdekis & Arimoto 
1999 1, but the behaviour of these diagnostics in stellar popu- 
lations with extended star formation histories is largely unex- 
plored. 

An alternative to using minor spectral features is to use the 
shape of the continuum over a wider wavelength range. The 
two degenerate stellar populations discussed in Sect. 4.2 (a 13 
Gyr old, r = 15 Gyr population and a 3 Gyr old, r = 1 Gyr 
population) have nearly identical continua throughout the opti- 
cal and near-IR part of the spectrum. As one moves to shorter 
wavelengths, this is however no longer the case. In the mid- to 
far-ultraviolet range (1000-3000 A), the slopes of the continua 
for these two populations differ significantly. Hence, the addi- 
tion of UV data could substantially improve the situation. Even 
the inclusion of U -band photometry would be useful. 

In this paper, EW(Ha) has admittedly been used in a very 
crude way, mainly due to the difficulties involved in sensibly 
weighting emission line data against photometry. A more in- 
telligent age fitting procedure could possibly improve the con- 
straints. The effects of underlying absorption and selective ex- 
tinction (Calzetti et al. 1994 1 would then however also have to 
be more carefully considered. 

Even if the ages turned out to be very low, it would be pre- 
mature to claim a contradiction with the underlying ACDM 
cosmology, since current simulations still have difficulties in 
forming realistic disk galaxies (e.g. Primack 2004 1. A partic- 
ular concern is that there may exist mechanisms which de- 
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lay star formation, so that while mass assembly of LSBGs 
takes places early in the history of the universe, in accord with 
ACDM, significant star formation does not occur until much 
later. Photoionization by the extragalactic UV background has 
been advocated as one way to suppress star formation in very 
low-mass galaxies (Dong et al. 2003 1, but would not be effi- 
cient for objects in the range of total masses relevant for our 
galaxies (probably ~ 10 s — 10 11 M & after inclusion of dark mat- 
ter). A more viable mechanism is perhaps suppressed star for- 
mation due to high angular momentum. One popular explana- 
tion for the low star formation efficiency in LSBG is that these 
galaxies formed inside dark matter halos with unusually high 
spin parameters (e.g. Dalcanton et al ll997l Jimenez et al. 119981 
Boissier et al. 2003 1. A similar scheme has also been suggested 
to explain the extremely low star formation rates in the hypo- 
thetical dark galaxy population (e.g. Jimenez et al. U997> . Until 
processes like these, which could substantially affect the SFH 
of stellar populations, have been considered in the predictions 
of the average formation redshifts for the stars in galaxies, the 
assumed cosmology cannot really be blamed for any potential 
age discrepancy of the type discussed here. 



6. Summary 

In a critical assessment of the possibility to derive properties of 
the bluest low surface brightness galaxies from optical/near-IR 
broadband photometry and Ha emission line data, we find that: 

- The global EW(Ha) observed in the bluest LSBGs are too 
low to be consistent with recently proposed scenarios in- 
volving constant or increasing star formation rates over cos- 
mological time scales (Boissier et al. 2003 ), unless the frac- 
tion of high-mass stars is significantly lower than predicted 
by the Salpeter IMF with M np = 120 M . A combination of 
EW(Ha) and chemical abundances may provide additional 
constraints on both the star formation history and IMF, and 
should be investigated further. 

- The ages of these objects are poorly constrained without 
prior knowledge of the star formation history. In fact, cur- 
rent observations cannot rule out the possibility that these 
objects formed as recently as 1-2 Gyr ago. Additional UV 
observations may however improve the age estimates. 

- If low ages are typical for the bluest LSBGs, this would 
contradict current models of hierarchical galaxy formation 
in the ACDM cosmology, since only a small number of 
galaxies with stellar masses typical for our objects are pre- 
dicted to form this late in the history of the universe. 

- Contrary to common methodology when estimating the 
masses of stellar populations, we note that the uncertain- 
ties associated with stellar M/L are not necessarily smaller 
in the near-IR filters. Instead, the /-band appears to be op- 
timal. 
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